+ 1 formalism in Einstein-Gauss-Bonnet gravity 
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Towards the investigation of the full dynamics in higher-dimensional and/or stringy gravitational 
model, we present the basic equations of the Einstein-Gauss-Bonnet gravity theory. We show (A^-|-l)- 
dimensional version of the ADM decomposition including Gauss-Bonnet terms, which shall be the 
standard approach to treat the space-time as a Cauchy problem. Due to the quasi-linear property 
of the Gauss-Bonnet gravity, we find that the evolution equations can be in a treatable form in 
numerics. We also show the conformally-transformed constraint equations for constructing an initial 
data. We discuss how the constraints can be simplified by tuning the powers of conformal factors. 
Our equations can be used both for timelike and spacelike foliations. 

PACS numbers: 04.20.Ex, 04.25. D-, 04.50.-h, 11.25.Wx 



I. INTRODUCTION 

General relativity (GR) has been tested with many ex- 
periments and observations both in the strong and weak 
gravitational field regimes (see e.g. and none of them 
are contradictory to GR. However, the theory also pre- 
dicts the appearance of the spacetime singularities under 
natural conditions [1, 0| , which also indicates that GR is 
still incomplete as a physics theory that describes whole 
of the gravity and the spacetime structure. 

We expect that the true fundamental theory will re- 
solve these theoretical problems. Up to now, several 
quantum theories of gravity have been proposed. Among 
them superstring/M-theory, formulated in higher dimen- 
sional spacetime, is the most promising candidate. In 
these prospect, a considerable number of studies con- 
cerned with gravitational phenomena and cosmology 
have been made in the string theoretical framework be- 
yond GR. 

Since the present knowledge is still far from under- 
standing the full aspects of the string theory, several 
kinds of approaches, which are fundamentally approx- 
imation, are usually taken. Among them perturbative 
approach plays important roles. There are two particular 
parameters which characterize the system in the super- 
string theory. One is string coupling parameter = e"^, 
where (j> is the dilaton field. The other is the inverse 
string tension a'. When the tension is strong (i.e., small 
a') compared to the energy scale of the system, it is diffi- 
cult to excite strings and the size of the strings becomes 
small enough to be regarded as particles in the zeroth 
order approximation. In this limit GR (with other light 
fields) is recovered. This is called Q;'-expansion[l|. 

In the higher order terms of a' , curvature corrections 
appear. The Gauss-Bonnet (GB) term is the next lead- 
ing order of the a'-expansion of type IIB superstring 
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theory 3, Hi, and has nice properties such that it is ghost- 
free combinations and does not give higher derivative 
equations but an ordinary set of equations with up to 
second derivative in spite of the higher curvature combi- 
nations. 

The models with the GB term and/or other higher cur- 
vature terms have been intensively studied in the high en- 
ergy physics. One of them is a series of studies in string 
cosmology. The pre-big bang scenario is a fantastic 
scenario which tries to avoid the Big-Bang singularity by 
making use of T-dualityQ (or scale factor duality). Fur- 
thermore, the pre-big bang scenario gives a natural in- 
flation mechanism, since the solution in the pre-big bang 
phase is inflating from the beginning at least in the string 
frame. Although these analyses show that the singularity 
problem has not been resolved yet completely, there are 
some cosmological solutions which do not start from an 
initial sing ularity[i,[i3,[ili. 

In addition to cosmology, the string effects can be also 
seen in the study of black hole physics. As the size of 
a black hole becomes small, and the curvature around 
the black hole becomes large; it is expected that the 
curvature corrections cannot be negligible. The singu- 
larity inside of the event horizon would be also modi- 
fied or even disappeared by string effects. For these rea- 
sons, the static or stationary black hole solutions in ef- 
fective string theories were investigated both in the sys- 
tems without higher curvature terms p^. ITsj and with 
such terms [3, Besides the static or stationary 
solutions, there are some dynamical solutions which are 
motivated by gravitational collapse [iTj. 

These analysis arc performed on the assumption of 
highly symmetric spacetime because the system is much 
more complicated than that in GR. To obtain deeper un- 
derstanding of the early stage of the universe, singularity, 
and/or black holes, we should consider less symmetric 
and/or dynamical spacetime; the analyses requires the 
direct numerical integration of the equations. None of 
the fully dynamical simulations in GB gravity has been 
performed. 

In this article, we present the basic equations of 
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the Einstein-GB gravity theory. We show {N + 1)- 
dimensional version of the ADM decomposition, which is 
the standard approach to treat the spacetime as a Cauchy 
problem. The topic was first discussed by Choquet- 
Bruhat [31 , but the full set of equations and the method- 
ology have not yet been presented. In 4-dimensional 
GR, numerical simulations of binary compact objects are 
available in the last years, and many groups apply the 
modified ADM equations in order to obtain long-term 
stable and accurate simulations. However, such mod- 
ifications are depend on the problem to consider, and 
the "robustest" formulation is not yet known (see e.g. 
p^). Therefore, as the first step, we in this paper just 
present the fundamental space-time decomposition of the 
GB equations, focusing on GB term. 

The ADM decomposition is supposed to construct the 
spacetime with foliations of the constant-time hypersur- 
faces. This method can be also applied to study the 
brane-world modeljlj, which states the visible space- 
time is embedded in higher dimensional "bulk" space- 
time. As was first investigated by Chamblin et al [2l|, 
it is possible to study the "bulk" structure by switching 
the normal vector of the hypersurface from timclike to 
spacelike. We therefore present all set of equations for 
both cases for future convenience. 

The outline of this paper is as follows. In Sec. HH we 
show that the set of equations are divided into two con- 
straints and evolution equations along to the standard 
procedure. In Sec. Illli we present the conformal ap- 
proach to solve the constraints which shall be used for 
preparing an initial data. In Sec. IIVI wc show the dy- 
namical equations in GR and in GB theory separately. 
Sec. |V] is devoted to the discussions and summary. We 
think these expressions arc useful for future dynamical 
investigations. 



II. (iV + l)-DECOMPOSITION IN 
EINSTEIN-GAUSS-BONNET GRAVITY 



A. Model and Basic Equations 

We start from the Einstein-Gauss-Bonnet action in 
(A'^ -|- l)-dimensional spacetime {M.,g^v) which is de- 
scribed as 1271: 
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(1) 



(2) 



where is the {N + l)-dimensional gravitational con- 
stant, TZ, TZf_,y, TZ^i,pa and £matter are the (A^ + 1)- 
dimensional scalar curvature, Ricci tensor, Riemann cur- 
vature and the matter Lagrangian, respectively This 



action reproduces the standard (A'^-l- l)-dimensional Ein- 
stein gravity, if we set the coupling constant acs (> 0) 
equals to zero. 

The action ([1]) gives the gravitational equation as 



where 



+'R^'^^'Rual3'y - -^gtj.uCGB, 



(3) 
(4) 

(5) 
(6) 



B. Projections to Hypersurface 

In order to investigate the space-time structure as the 
foliations of the A^-dimensional (spacelikc or timclike) hy- 
persurface E, we introduce the projection operator to E 
as 

-Lfiu = gtiv - enf^n^, (7) 

where is the unit-normal vector to S with n^rt'' = 
e, with which we define is timelike (if e = —1) or 
spacelikc (if e = 1). Therefore, E is spacelikc (timclike) 
if is timclike (spacelike). 

The projections of the gravitational equation ([3]) give 
the following three equations: 

{Q^ly + acBTi-i^iu) n" T^y n" = K^p, (8) 



(9) 

{QfLV + OLQB'Hfj.v) -L'^ -L'^o- = T^LV -L'^ -L'^o- = I^^Spa, 

(10) 

where we defined the components of the energy- 
momentum tensor as 

T^iy ^ pn^n^ + Jf^n^ + J„n^ + S^„, (11) 

and wc also define T = ep + 5""^ for later convenience. 

Projection of the (A^-|- l)-dimensional Riemann tensor 
onto the A'-dimensional hypersurface can be written as 

R^upa ^(A^pAT^yjj K^fjKup')^ (12) 

n^p^s 1-% r'p n' = -2i?[^AV]p, (13) 

7^„/37^ nf' = f^K^p + if^^ if"^ , (14) 

where Rp,vpa is the Riemann tensor of the induced met- 
ric 7^,y(= -L^i/), Dp is the covariant differentiation with 
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respect to 7^1/, £n denotes the Lie derivative in the n- 
dircction, and K^^^ is the extrinsic curvature defined as 

K^,, = -]^£,,-t^, = -l.''^l.P^V^np (15) 
Eq. (fT^ is called the Gauss equation, and the contraction 



of ((T3l) of ji and v gives the Codacci equation. 

Using these projections, the {N + l)-dimensional Rie- 
mann curvature and its contractions (the Ricci tensor 
and scalar curvature) are described by the iV-dimensional 
variables on the hypersurface S as@ 



^p'^^al^ua^ p 



^u^pl^pa^ (J ^u^al^ pa^ p ^^p^p^n^ua ^p^^a £ i^p 
l^vl^p^ per ~t~ f^u^a £ ppt 



(16) 



KKp, - 2Kp^K% + Up {Da^Kl - D,K) + {D^K"^ - DpK) 



K^pK''^ + e£nKp 



(17) 
(18) 



where K — K"^. 

Substituting these relations into the field equation ([3]) or ([5|)- ([TU)) . we find the equations are decomposed as (a) the 
Hamiltonian constraint equation: 

M + UGB (M^ - UlabM"''' + MabcdM''''^'') = -2eK^pH + 2A , 



(b) the momentum constraint equation: 

and (c) the evolution equations for 7^ : 

Af,, - iM7,, - e{~K,aK-^ + j,,KabK^' - £„K,^ + j.^^^" £,,Kab) 



2a 



GB 



H,, + e{M£nK,j - 2M^''£nKaj - 2A//i:„ifa» - W,j'"'£nKab] 



respectively, where 

Mijki = Rijki — £{KikKji ~ KuKjk), 
Nhj = 7"'M„jfc = R^J - e{KK^j - K^aK"^), 

M = r'Mab - i? - £(A-2 - KabK''''), 
Nijk = DiKjk — DjKik, 

N^ = Y'''Na^b - DaK/' " D,K, 

H,j = MAU, - 2{A'haM"-^ + Ar'Ahajb) + M.^bcM^'^'' 



-2e 



kl 



-^7y [M2 - 4MafcM''^ + A'labcdM'^'''^] 

-£7„- [KabK''''M ~ 2AdabK'"'K^ - 2NaN'' + NabcN""'"] , 

M7„-7'=' - 2M„-7'=' - 27,jA/'=' + 2A.f,,jb7'''=7''- 



(19) 
(20) 



(21) 

(22) 
(23) 
(24) 
(25) 
(26) 



(27) 
(28) 
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Wc remark that the terms of £nKij appear only in the 
hnear form in (|2ip . This is due to the quasi-hnear prop- 
erty of the GB gravity. 

C. Cauchy approach 

Using the Bianchi identity, Choquet-Bruhat [l8l | 
showed that the set of equations forms the first-class sys- 
tem as the same as in GR, that is a space-like hypersur- 
face which satisfies the constraints (|T9|) and ([20|) will also 
satisfies the constraints after the evolution using ([2T|) . 

The most standard procedures for following the dy- 
namics of space-time consist from the three steps: 
(i) solve the constraint equations (|19|) and (pUjl for 
{"^ij , Kij , Ji) on = 0) and prepare them as the 
initial data, (ii) evolve {jij, Kij, pn, Ji) using (|21[) and 
the matter equations, and (iii) monitor the accuracy of 
the evolutions by checking constraint equations on the 
evolved 

In the case of seeking the "dynamics" along a space-like 
direction x such as a study of the "bulk" structure in the 
brane-world model, the above strategy can be switched 
to the evolution in x-coordinate instead of t (using the 
set of equations of e = +1). The initial data, in this 
case, is a time-like hypersurface which should satisfy the 
constraints. Such an initial data can be obtained either 
by solving the dynamics of "brane" part or by taking 
double Wick-rotation after the above step (i) , depending 
on the models and motivations. 

In the following sections, we describe a way of solv- 
ing constraints fSec. Illlj ). and a way of solving evolution 
equations fSec. lIVp. 

III. CONFORMAL APPROACH TO SOLVE THE 
CONSTRAINTS 

A. "Conformal Approach" 

In order to prepare an initial data for dynamical evo- 
lution, we have to solve two constraints, (fT9|) and (f20|) . 
The standard approach is to apply a conformal transfor- 
mation on the initial hypersurface [1^. The idea is that 
introducing a conformal factor ip between the initial trial 
metric jij and the solution 7^ , as 

Ty—V''"^-, 7'^ =^"""f^ (29) 

where to is a constant, and solve for ip so as to the solu- 
tion satisfies the constraints. 

For A^-dimensional spacetime, Ricci scalar is trans- 



formed as 

+ {N -1)[2- {N -2)m]mil)~'^{bi)fy (30) 

- (iV - 2)mip-^D,D.jiP 

+ {N - 2)to(to -I- l)%p~'^b,ipDj^ 

-TO[(iV-2)m- l]^-2(^V)'7y, (31) 

+ mij-^u [bjbki' - (to + l)^~^b,i^Dk^] 

- mij-^^,k [bjbii' - {m + i)^-^bji^biij] 

+ mi'-^^jk [DMi' - {m + 1)^-1 AV'A'/'] 

- miP^^-ji [b,bki^ - (to + l)ir^b^i^bkip\ 

+ m'^~^{b^ f{%ajk - %klji)]- (32) 

Regarding to the extrinsic curvature, we decompose 
Kij into its trace part, K = Y^Kij, and the traceless 
part, Aij = K jj — j^jijK, and assume the conformal 
transformation [28| as 

A, = yj'A,,, A'^ =i;'-^"-A'^, (33) 
K = iP^'K, (34) 

where £ and r are constants. For the matter terms, we 
also assume the relations p = ^~Pp and J* = tjj~'^J'^, 
where p and q are constants, while we regard the cosmo- 
logical constant is common to the both flames, A = A. 

Up to here, the powers of conformal transformation, 
£,m,T,p and q are not yet specified. Note that in the 
standard 3-dimensional initial-data construction cases, 
the combination of m = 2, £ = —2, r = 0, p = 5 and 
q = 10 is preferred since this simplifies the equations. 
We also remark that if we chose t = £ — 2to, then the 
extrinsic curvature can be transformed as Kij = ^^Kij 
and K'^ = V'^"'*™i'''^ . 

B. Hamiltonian constraint 

Using these equations, the Hamiltonian constraint 
equation turns to be 
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2{N - l)mbaD''i^ -{N -1)[2-{N - 2)m] m^il^fi'^ 

= R^P- ^X^£^2m+2r+i^.2 ^ e^,-2m+2f +i^^^^ab ^ 2eK^ p^^'P - 2A + aGse^A^^+i. (35) 

We will show the explicit form of the GB part Q = - AMabM'^'' + MabcdM"'""^ in Appendix El At this moment, 
we observe that (|35|l can be simplified in the following two ways. 

(A) If we specify t ~ l~ 2m and m = 2/{N — 2), then (I35p becomes 

^i^^D.D-i; = e4'^'+'-^/<-^-'\K^ - KabK'") + 2eK'p^-P -2A + aaBQi^^^^'^''-^^ (36) 

In the case of the Einstein gravity {acB = 0) with A = 0, the combination £ = 2/{N — 2) and p = —1 makes 
the RHS of ([55)) linear. If we choose I — —2, which will make the momentum constraint simpler as we see later, 
pop also remains as a simple equation. 

(B) If we specify t — and m = 2/ {N — 2), then ([55)) becomes 

(37) 



C. Momentum constraint 

In order to express the momentum constraint equation in a tractable form, additionally to the variables in Sec. lIII K\ 
we introduce the transverse traceless part and the longitudinal part of A^^ as 

DjAi^rj. = 0, (38) 
A'l = A"^ - i^-r, (39) 

respectively. Since the conformal transformation of the divergence DjA^-' becomes 

DjA^^ ^ 7/.^-2'»{z)ji/' + V^-i [i + m{N - 2)]i/'i)jV;}, (40) 
the momentum constraint, ([20[) . can be written as 

ij'-'^"'DaA.^l + [£ + (TV - 2)m] ~ ^X_Ij)^(^-AO + 2aGBS, = (41) 

I 



We will show the explicit form of the GB part in 
Appendix jB) meanwhile we proceed the standard proce- 
dure. That is, introducing a vector potential for 
as 

A'l = D'W^ + - —fWkW''. (42) 



Since the divergence of becomes 
» N — 2 

DjA'l = D^i:*"!^' + ^^D'L'fcW^'' + R\W'', (43) 
the momentum constraint, ()4ip . becomes 



+ ^^^AT^feVK'^ + + V"' + (iV - 2)m] (^"I¥^ + ^^W^'^ - -^r''DkW^)%^Da^ 

Similarly to the case of the Hamiltonian constraint equation, we consider the two cases. 
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(A) If we specify t = I — 2m and m = 2/ (TV — 2), then (|44p becomes 

^ N-2 
iV- 1 



N 



£- 



N-2 



{D,ij)K + D,K 



(45) 



In the case of the Einstein gravity {acB = 0), the choice of ^ = — 2 cancels the mixing term between ip and 
W*. Further, when acB = 0, we have a chance to make two constraint equations, and decouple by 
assuming K = and q = 8/{N — 2) + 2. However, when aas 0, this decoupling feature is no longer available, 
since the term includes '0-related terms as we see in eq. (jB7l) in Appendix [BJ 



(B) If we specify t = and m = 2/{N -2), then (gH) becomes 

A^- 2 



N-2)- 



.N- 1 
N 



(46) 



For the Einstein gravity {acB = 0), the choice of = —2 again cancels the mixing term between ij] and W. 
The decoupling feature between ([37|) and (f46|) is available when acs = 0, =const. and q = 8/{N — 2) + 2. 
However, when acB 7^ 0, this decoupling feature is no longer available, since the term includes '(/'-related 
terms as we see in eq. (jB13p in Appendix |B] 



D. Procedures 

For the readers' convenience, we summarize the above 
procedure briefly. The initial data, {jij , Kij , p, J^), can 
be constructed by solving the Hamiltonian constraint, 
(fT9|) . and the momentum constraint equations, ([20| . This 
can be done by the following steps. 

1. Give the initial assumption (trial values) for 7^ , K, 
ATT^ and p, J. 



we obtain the solution jij , 
the constraints. 



Kij , p, J* , which satisfy 



2. Solve ([35]) and (|44| for -0 and by fixing the 
exponent £, m, T,p and q. 

3. By the following inverse conformal transformations. 



(47) 

] 

(48) 

(49) 
(50) 



E. Momentarily static case 



The most easiest construction of an initial data may 
be under the assumption of the momentarily static (or 
time-symmetric) situation, Kij = Ji = 0. In such a 
case, the momentum constraint becomes trivial, and the 
Hamiltonian constraint (1351) can be reduced as 



2{N -l)mf)ab''^- (A - 1)[2- {N ~2)m]m{bijfi;'^ = R^j + 2eK^ pi^'^ ~ 2A aGsOV-^^+S (51) 
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where 



R~{N -2) [{N - 3)m - 2] mi^^^iD^jf DaD°"4^ 



R"'' + iN - 2)m(m + l)!/^"^^"^!)^ 
+ {N- 1)2TO^ [{N - 4)™ - 4] il)-^{bi}}f - 2^}-^ [{N ~ 4)to - 2] {DijjfR 
- 8(m + l)^-^R''''ba^bbi^^ + ^-^""Rgb, 

I 



(52) 



where 
{D-n\ 



{D-n)\ 
{D-m - 1)! 



= {D~n){D-n-l) ■ ■ ■ (D-m), 

(53) 



'ihcd 



and RgB — R — iRabR"' + RabcdR"""^ 

When = 3, 8 smiply becomes <d = ijj^'*'"''RGB so 
that ([5T|) will be reduced as 

Sbab"^ ^ R^ + 2eK^p^p-P -2A + aGBRcB^J'^, (54) 

for the choice of m = 2. 

However, in general N, it is hard to find out an ap- 
propriate m which simplifies the equation, even if jij is 
taken to be the flat space-time that reduces Q as 

e = {N- 3)™^^"''""^ 

X |4(7v - 2) [{babyjf - {babbtij){b''b''^j) 

+ 4:{N - 2)[{N - 3)m ~ 2]^l;-'^{bijfbab''^ 
+ 8{N - 2)(m+ l)^-i(£)»(£)''V)£>a£'bV' 
(iV - f )2m [(iV - 4)m - 4] i^-^ (bt/j)-^ I . (55) 



Roughly speaking, Eq. ((55|) is a quadratic equation 
with respect to the second order derivative of ip, which 
means that there are two roots in general when a set of 
trial values ■jij and p on the hypersurface are given. The 
meaning of the existence of two solutions is more clearly 
understood by assuming p — and q = in the confor- 
mal transformation of the matter field. In this case, two 
different 7^ are obtained through the different conformal 
transformations even for the same matter field distribu- 
tions. 

The non-uniqueness of the solution is due to the higher 
curvature combination of the GB terms. The existence 
of two solutions can be also seen for the black hole solu- 
tions in the Einstein-GB theory [11 [11]. In the black hole 
case, the gravitational equations reduce to a quadratic 
equation under the suitable choice of a metric function. 
The two roots of the equation correspond to the differ- 
ent black hole solutions. In the acB limit, one 
of the black hole solutions approaches to the solutions 



in GR while the metric component of the other solution 
diverges and there is no counterpart in GR. Hence the 
former is classified into the GR branch and the latter 
into the GB branch. In some mass parameter regions of 
the black hole, the metric components of these branches 
coincide at the certain radius. It is a multiple root of 
the field equation. Solving the field equations beyond 
this radius, we find that the metric becomes imaginary 
and this region is physically irrelevant. The spacetime 
becomes singular at this radius and it is called branch 
singularity, which is a new ingredients of GB gravity. 

When we solve the Hamiltonian constraint equation 
(|55p . the similar situation may occur. The conformal 
factor tp becomes imaginary and the solution "fij is also 
imaginary in some regions on the hypersurface. The 
boundary of this region corresponds to the singularity. 
The situation is same also in the non-time symmetric 
case. 



IV. DYNAMICAL EQUATIONS 
A. Dynamical equations in General Relativity 

Equations in GR are obtained by putting acB = 0. 
The evolution equations of (A^-f l)-dimensional ADM for- 
mulation are presented in [24| . but here we generalize 
them to the set of equations both for spacelike and time- 
like foliations. 

Firstly we rewrite the dynamical equations (|2ip by in- 
troducing the metric components as 

ds^ = ea^{dx°f + -f^J{dx' + f3'dx"){dx^ + (3^dx°) 
= {ea^ + (3a(i''){dx°f + 2/3^dx'dx° + j,jdx'dx^, 



(56) 



where a and /3* are the lapse and shift functions, respec- 
tively. The components of the normal vector, then, are 

= (ea, 0, ••• , 0), = -(1, -/3'). (57) 

a 

In the matrix from, the {N + 1)- and N- metrics are 
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expressed as 



1 -p' 



and 



7'^" = 





7*^' 



(58) 



(59) 



(60) 



Trace of Eq. (HH) minus half of Eq. ^ gives 

^2 



iV- 1 



r + eA. (61) 



This is equivalent to the trace of Einstein equation ([3]) 
with the projection p^ . Substituting Eq. (pT|) into 
Eq. dUl), we find 



/ K 



en^ ( 5',:-,- - 



1 \ 2£ 



The extrinsic curvature and its Lie derivative are ex- 
pressed as 

K,, = {-^oJ^J + D,p, + D,Pj) , (63) 



a 



Ka.D.p'' - Ka.D.P'') 



(64) 



respectively. 

With the metric components ([SG]) . the equations ([63]) 
and (16211 are 



-2aify + A/3j, (65) 



doK .,j - 



+ iD,P'')K,a + (Ar )i^a 



1 



N - 1 



2as 
N - 1 

(66) 



If we have matter, we need to evolve them together 
with metric. The dynamical equations for matter terms 
can be derived from the conservation equation, V^T^^ = 
0. 



B. Dynamical equations in Gauss-Bonnet gravity 

With the Gauss-Bonnet terms, the evolution equation, 
PT|) . cannot be expressed explicitly for each £nKij . That 
is, (PT|) is rewritten as 



(1 + 2aGBM)£nK^j - (7^7''^ + 2aGBW,/'')£nKab - 8aGsM(/£„i^|„|j) 



:M7,. 



(67) 



and the second and third terms in r.h.s include the for the Gauss-Bonnet gravity. We expect that in the most 
linearly-coupled terms between £nKij. Therefore, in 
an actual simulation, we have to extract each evolution 
equation of Kij using a matrix form of Eq. (I67p like 



k = Ak + b 



(68) 



where k = {£nKii, £nKi2, • • • , £nKfqMY' ^'^'^ t) are 
appropriate matrix and vector derived from Eq. (|67|) . 

The procedure of the inverting the matrix (1 — ^) is 
technically available, but the invertibility of the matrix 
is not generally guaranteed at this moment. In the case 
of the standard ADM foliation in 4-dimensional Einstein 
equations, the continuity of the time evolutions depends 
on the models and the choice of gauge conditions for the 
lapse function and shift vectors. If the combination is not 
appropriate, then the foliation hits the singularity which 
stops the evolution. The similar obstacle may exist also 



cases Eq. (j68p is invertible for Kij but we cannot deny 
the pathological cases which depend on the models and 
gauge conditions. Such a study must be done together 
with actual numerical integrations in the future. 



We obtain, however, the similar form of equation by 
introducing the {N + l)-dimensional Weyl curvature [25l 
[26j . It is useful for the discussion of the brane- world 
cosmology although it is not written by only the A'^- 
dimensional quantities, and we can not adopt it as the 
evolution equation directly. Trace of Eq. (|2ip minus half 
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of Eq. dUl) gives 
(N-l) 



M + KabK'^' + h''''£nKab 



+2(n - 3)aaB 



+MKabK''^ - 2K)K{m\ - 2NaN'' 

+NabcN''''^ + Mh^^lnKab - 2Ar^£nKab 

= -£K^T + eA 



(69) 



By the last term in the l.h.s of this equation, the term 
of the Lie derivative cannot be expressed by other term 
exphcitly. 

Let us rewrite the dynamical equation (|2ip in a dif- 
ferent form. From Eqs. IHl), dTT]) and (0 with the 
decomposition of the Riemann tensor as 



V — 



2 

- jV(j\r_ ^^ 9p^[p9a\u'Tl + C^^pCT , (70) 



where Cp^pa is the (N+l)-dimensional Weyl curvature, 
we find 



iV - 1 
£nKij — £ ^ _ 2 



N -2 



M, 



(71) 



where 



E,, := Cp^p, n^^nP^^j'^j . (72) 



However, because Eq. ([7T]) is a trace free equation 
£nKpi, cannot be fixed by Eq. ([7T|) . Inserting Eq. ([7T| 
into Eq. we find 



h°'''£nKab=-^M-KabK'''' 



-^('^^T-A) + M^^^il^/,(73) 



where 



.r . 2(7V-2)(iV-3) 



(74) 



N -6 
4:{N-2)' 



■M 



#2 



N 

JL^M^bM'^' + ^MabcdM-'"" 



NaN'' + iNabcN' 



2(A^~1) 
iV- 2 



MabE-' 



From Eq. ^ with Eq. jTSl), we then find 



N - 1 



N -2 



K — 

" J NU 

2e{N - 2,)aGB 
NU 



{n^T - A)7y 



^7, 



(76) 



Since Eij is (7V-|-l)-dimensional quantity, £nKij can not 
be evaluated by the valuables on A''-dimensional hyper- 
surface with this equation. This means that the Eq. ([75]) 
can not be used for the full dynamics as we mentioned. 
For example, however, for the Friedmann brane-world 
model, where the constant time slice of the timelike hy- 
persurface is homogeneous and isotropic, Eij can be writ- 
ten using quantities on the hypersurface. For such lim- 
ited cases where the term Eij can be evaluated on the 
hypersurface, the Eq. ([75]) is useful for simplifyng the 
situation. 



V. DISCUSSION 

With the aim of numerical investigations of spacetime 
dynamics in higher-dimensional and/or higher-curvature 
gravity models, we presented the basic equations of the 
Einstein-Gauss-Bonnet gravity theory. 

We show the {N -f l)-dimensional decomposition of 
the basic equations, in order to treat the space-time as a 
Cauchy problem. With the aim of investigations of bulk 
spacetime in recent brane-world models, we also prepared 
the equations both for timelike and spacelike foliations. 
The equations can be separated to the constraints (the 
Hamiltonian constraint and the momentum constraint) 
and the evolution equations. 

Two constraints should be solved for constructing an 
initial data. By showing the conformally-transformed 
constraint equations, we discussed how the constraints 
can be simplified by tuning the powers of conformal fac- 
tors. If we have Gauss-Bonnet terms, however, the equa- 
tions still remain in a complicated style. 

For the evolution equations, we find that £nKij com- 
ponents are coupled. However, this mixture is only up to 
the linear order due to the quasi-linear property of the 
Gauss-Bonnet terms, so that the equations can be in a 
treatable form in numerics. 

We are now developing our numerical code and hope 
to present some results elsewhere near future. 



APPENDIX A: GB PART OF HAMILTONIAN 
CONSTRAINT EQUATION 



In Eq. f35)) , the GB part of the Hamiltonian constraint 



equation 9 = A/^ - AMabM"^ + MabcdM"''""' is not writ- 
(75) ten explicitly. It becomes as 
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e = (iV - 3)to?/'"'""<^ 4(7V - 2)mTp-^ (DaD^tl^y - DaDbiljD" D^'i}^ 



Ai)-^ M -{N - 2) [{N - 3)m - 2] m^j-^ Da^l^D^i^ 



+ {N- \)2Tr? {{N - 4)m - 4] ^-^{ba^jt>''i)f - 2^^"^ [(A^ - 4)to - 2] MD^jD'^i' 



where 



T = R-e 



N ~1 



^ ij ^ 



^ijkl — Rijkl 



N 
N - 1 



When T = £ — 2m and m = 2/{N — 2), which corresponds to the case (A) [Eq. ([M]) ] 



2(7V-3) 2 
7V-2 ^ ^ ^ 



47/-"^ (M + 2tir^Daijb''^)bab''tp + Si/-"^ (iir'' + 



2iV 
iV-2 



8N{N - 1) 



(iV-2)2 " ^ ■ iV-2 



8N 
N -2 



where 



(Al) 



(A2) 
(A3) 

(A4) 



T = R- £V2^-4/(^-2)(i^2 _ KabK'^"), 

t,, = R,, - e^^'-''/^^-^\KK,, - k,aK^^), 

When r = 0, m = 2/(A^ — 2), which corresponds to the case (B) [Eq. ([57]) ]. Q is expressed as Eq. (|A5|) and 
'N -1 



(A5) 



(A6) 
(A7) 
(A8) 



T = R-e 



Tij — Rij 



' ijkl ^ijkl 



N 

'N -1 



^4/(iV-2)^2 _ ^2£-4/(iV-2)^^^^a 



(A9) 
(AlO) 

(All) 
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APPENDIX B: GB PART OF MOMENTUM CONSTRAINT EQUATION 

In Eq. (j44p . the GB part of the HamUtonian constraint equation is not written exphcitly. It becomes as 
Ei = V-^-^™!^ - 2{N - 3)mip-^DbD''iP ~ {N ~ 3)m[{N - 4)m + 2]^|J-'^ DbiPD''^P 



iV2 

_^_^^-4m|2^ah _ 2(iV - 3)TOV'-ii)^i)°i/' " 2(iV - l)m(m + b^ipb^ip 



- ^^^^^^ [(TV - 4)m + 2f + r] eV-^'-^^+^-^xAV'i'ai", 

2(jV^ 3) g^2£-4m+r j,^ j^^b^^a _ 2 [( - 6)m + 3^] eV'^^-^"-! AV-i^bi^i" , (Bl) 



where 



7^, = |[(A^-3)m + ^]V'^-4™-ii,°AV'- ^^^^V'"^"+^(A-ft^ + TV'"^-^AV')|^ 

^(^^_^^^-2m+r-l^^^^ + ^i^X_J^^-2™+r^^^ _ 2 [(TV - 3)m + ^^^-^^-l^^b^^^j^^a 

-2(m - l)i,'-^-^-\Kbb,i, - A,bba^)it' + 2(m - Av^ifcci?,™', (B2) 

2?. - ^ ~ ^ ^\ n^-^-+"-^(Ai^ + Tir'kb,^) 



~ 2m [(iV^ „ 6iV + 7)m+ (iV - 3)^] V/"**"'"^ AV'^^ j A-D'^V' 

-| ^^^^ ^^^ m^-^'"+-^(Aj^ + Ti^-'kbai^) 

~ 2m[(iV2 _ + 5)„j + (AT _ 2)(^ - 2)] V^^^-^^-^^j^^^b^j/jiT)^^ 
+2(iV - 3)m(m - £)^^-4m-2(^^^^^^ _ i,, AV')^'^"'/', (B3) 

^(1) = 2|-;^^^^m(TO + l)V'-2'"+2^-2p^^ _^ Tljj-^kba^j){b''k + T^J-^KWi^) 

(N — 2)'^ ^ ^ ^ N — 2 ... 

+ ^ ^ ' m{m + l)?/>-2"+^-2£,^^^a^ ^ __„j[(7V2 - 4Ar + 5)m + 2] TV'~^"'+''"^ii:L'aV'^''V' 

- (m - ^)£V^'^"^'""^ia^ih'ic''V (B4) 



{N - 3)£V' 



2£— 4m+T^ byi"'' 



^^(3) ^ -m[(iV- 2)2(Ar_ 5)^2 _^ (TV - 2)3(^- 2)™+ (A^- 1)(3£- 2)]V'^-'''"-3i)„V-D''V' 
- 4^ [(^ - - 8)m + (iV2 - TV + 2)^]eV'^-2m+2r-i^2 

+ [(iV - 6)m + 3^] 

When T = £ — 2m and m = 2/(A^ — 2), which corresponds to the case (A) [Eq. (|l5)) ]. 



' ^ 1 Ar-2 ^ ^ N-2 



- 3iV + 4 



, 2(iV-3) 

^^£-8/(iV-2) L^ab _ 4(iV~3) 1^,^ _ 4jV(jV- l) 



.^2^-4/(Ar-2)^^a6 ^ 2e^2£-4/(iV-2)^a^^cfc I (£) _ ^^1^^) 



2(7V-3) 

+ 2e^3f-12/(JV-2)^^a^;,c(^^^^^ _ ^^^^^^ ^ ^ ^ ^(1) + A^^'^D.^K + ^(3) ^^^a 

4 



2(iV - 2)3 3^_12/(jV_2^ ^ 
7V2 



+ 



2(N-3) 
N 



2{N -6) 

N-2 



N-2 



iV 



2{N - 6) 



TV- 2 



where 

7^, = ^^-8/(^^-2) 



2(iV- 3) 



7V-2 



N-3 
N 



D 



k+{t 



N-2 



N-2 



N 



2{N -i) 
N-2 



R 



^-8/(W-2)-l 



DalpAbcR 



cab 

t 1 



^ \ N{N-2) ^ 



N-2 



{N - 3)e - 



2{N^-6N + 7) 



N-2 



^^8/(A^-2)f 4(iV-3) 1 

\ N ^ 



DaK 



N-2 



1 + 



6 



4riV- 3) 



{N - 2)2 
2 



^-^Dt^bA^ \D-b,^ 



/ _ 2 \ 
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iV- 2 



£-4/(7V-2)-2 



iV-2 



iV- 2 



1 



7V-2 



a6 



(BIO) 



^(2) ^ l^£-8/(iV-2) 1 _ 4(iV - 3) ^^^2fy^^fja 



N 



N -2 



(Bll) 



1 

3f 



iV- 2 
(A^2 _ TV + 2)^ + 
2(iV- 6) 



2(7V- 1) 



N ~2 



{N - 2){N^ - 8) 

^^2«-4/(Ar-2)-l^^^^ab 



£^2^-4-1^2 



(B12) 



When T = and m — 2/ {N — 2), which corresponds to the case (B) [Eq. (I46p ]. 



N -2 



(N-2Y 



]V2 ' 



.^4/(7V-2)^2 _ e^2^-4/(W-2)^^^^fcc ^^^a 



+ 



2(^-3)^,,. 



N 



_2(iV_2)3^^,_4/(A^-2)^^^^^. + ^^^^f^+ ^l^^'^-'/^^-'^-'^^^M^i'^. 
^ 2(A^^- 3) g^^2£-8/(iV^2)^^j^^fc^^a _ 3 



N V iV-2 
2(7V- 6) 



3^ + 



N-2 



(B13) 



where 



2(iV-3) 



iV- 2 



^/-8/(W-2)-l^^a^^^ 



iV-3 



-2\ 



TV 

2(iV-3) 



^-4/(iV-2)^^^ U 



Ar-2 



+2(£ 



7V-2 



«-8/(Ar-2) 



(B14) 



14 



V,, = 



2 {N'^-m + n 



N -2 



N 



{N-3)i 



^_4/(JV-2)-l^^^ 

2{N^-6N + 7) 



N~2 



N -2 



£- 



6 



(B15) 



= ^V^-4/(Ar-2)-2^^_^^a_^^4^-4/(Ar-2)-2^^^.^a^ 



7V-2 



8iV(iV - % i-8/(N-2)^3f)a^^f)b^J^ _ 2(jV - 3) / ^_V^2^-8/(^-2)-lj^i ia6 
(iV-2)2 ^ V aft AT V AT -2/ 



(B16) 



iV- 2 
1 



3f 



{N^ - N + 2)1 + 
2(iV-6) 



TV- 2 



2{N - 1)(7V2 - 8) 
7V-2 



,^£-4/(W-2)-1^2 



(BIS) 
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